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Abstract 

Quasi-elastic (z > 0.95) photo-production of ip{2S) mesons has been observed at 
HERA for photon-proton centre-of-mass energies in the range 40 to 160 GeV. The 
ip(2S) mesons were identified through their decays to £ + l~ and to J /ip tt + tt~ , where 
the J/ip subsequently decays to £ + £~, the lepton i being either a muon or an 
electron. The cross-section for quasi-elastic photoproduction was measured to be 
(18. 0± 2.8(stat) ± 3.0(syst)) nb at a photon-proton centre-of-mass energy of 80 GeV. 
The ratio of the *p(2S) to J ftp quasi-elastic cross-sections is (0.150 ± 0.027(stat) ± 
0.022(syst)). 
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1 Introduction 



At HERA the production of vector mesons is studied at photon-proton centre-of-mass energies 
(W 7 p) extending well beyond those reached in fixed target experiments. In particular, the 
production of J/ip mesons has been reported upon by HI |, g, | and ZEUS g, §, §. In this 
letter a study of the quasi-elastic photo-production of ip(2S), i.e. -0(3685), mesons in the photon- 
proton centre-of-mass energy range 40 < W 7P < 160 GeV is reported. Quasi-elastic is defined 
here as the kinematic region z > 0.95, where z = E^^s)/ i n the rest frame of the proton. 
The tp(2S) mesons are identified via the decays tp(2S) — > J/ip(— > £ + £~)it + it~ or ip(2S) — > £ + l~ , 
where I denotes a muon or an electron. 




Figure 1: Elastic production of charmonium 



At the collision energies studied here, the process 7p — > Vp (V = J ftp or ip(2S)) can be 
considered in three parts (see figure |l]): the formation of a cc quark pair from the photon, the 
interaction of the quark pair with the proton, and subsequently the formation of the charmonium 
state. In a QCD based model @, § of charmonium production described in section §, this is 
justified by arguing that, viewed in the proton rest frame, the cc quark pair creation from the 
photon and the formation of the charmonium state from the cc pair both take place over distances 
considerably greater than the proton radius. The cc quark pair is regarded as a colour dipolc 
which interacts with the proton. A prediction of the ratio of the cross-sections for ip(2S) and 
J/ifj mesons arising from this model is compared to the measured ratio, presented in section [?]. 



2 Kinematics 



The relevant kinematic quantities are defined in figure ||. For this study, it is demanded that 
neither the system Y nor the scattered positron e' are observed (see section ^|) , escaping in the 
+z (incoming proton) and — z (incoming positron) directions respectively. The system Y is the 
scattered proton in the elastic case, and otherwise can be considered as the dissociated proton. 

The photon-proton centre-of-mass energy, W 7P , is determined using the relation: 

'P 2E e [ ' 

where E e is the incoming positron energy, and the sum is over the energies, E, and the lon- 
gitudinal momenta, p z , of the decay products of the ip(2S) meson. The centre-of-mass energy 
of the positron-proton collision is denoted s. Neither the elasticity of the interaction, z, nor 
the photon virtuality, Q 2 , can be measured. However, the selection criteria restrict the accep- 
tance (see section |5|) approximately to the kinematic region z > 0.95, called here quasi-elastic, 
and Q 2 < 4 GeV 2 ; the median Q 2 is expected to be approximately 10 -4 GeV 2 . Cross-sections 
are evaluated and presented in the photo-production limit, Q 2 = GeV 2 , in section || All 
cross-sections are integrated over the momentum transfer variable t. 
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s = (P + k) 2 = 90200 GeV 2 

W 2 p = {q + P) 2 

Q 2 = -q 2 = —(k - k') 2 

y = P-q/P-k 

z = (-P/q- P 

t=(P- P') 2 



Figure 2: The kinematic variables describing ip(2S) production at HERA. The incoming positron 
and proton are denoted e and p respectively, the exchanged photon by 7, the outgoing positron 
by e', and the ip(2S) meson produced by ip(2S). The remainder of the final state is denoted Y. 
The associated four-momenta are indicated in brackets in the diagram. 

3 Experimental Conditions 

The data presented here were collected in 1994 and 1995 with the HI detector at HERA, when 
27.5 GeV positrons were in collision with 820 GeV protons. They correspond to an integrated 
luminosity of 2.7 pb" 1 in 1994 and 3.6 pb" 1 in 1995. 

The HI detector is described in detail elsewhere || and only components relevant to the 
present analysis are described here. Polar angles are defined with respect to the direction of the 
incoming proton beam. The Central Tracker (CT) surrounds the interaction region, covering 
the polar angle range 20° to 165°. Multiwire proportional chambers (CIP and COP) provide 
additional information for triggering in this region. Forward of the CT is the Forward Tracker 
(FT) covering the polar angular range 7° to 25°. Surrounding the trackers is a highly segmented 



liquid argon (LAr) sampling calorimeter [10] consisting of an inner electromagnetic section and 
an outer hadronic section. Polar angle coverage is from 4° to 153°. During the 1994 data taking 
a lead scintillator electromagnetic calorimeter (BEMC) covered the angular range 155° to 176° 
and was able to detect scattered positrons down to a Q 2 of 4 GeV 2 [pi]] . This calorimeter 
was replaced for 1995 data taking by a scintillating fibre and lead calorimeter (SpaCal) with 
coverage of the angular range 153° to 177.5°, thus extending the acceptance down to a Q 2 of 
1.2 GeV 2 [12]. The tracking chambers are immersed in a uniform axial magnetic field of 1.15 T 



provided by a superconducting magnet that surrounds the LAr calorimeter. Outside this lies 
the iron return yoke of the magnet which is instrumented to detect muons in the range 4° to 
171°. Luminosity is determined by measuring the rate of the reaction ep — > ejry. 

A first level trigger accepts events in which a i^(2S) meson is produced if a track of at least 
600 MeV momentum pointing towards an energy cluster greater than about 1 GeV in the LAr 
calorimeter is identified, or two tracks back-to-back in azimuth in the multiwire proportional 
chambers (CIP and COP) are found, or if there exists a muon signature in the instrumented 
iron of the return yoke. 



4 Identification of the i^{2S) 



The ip(2S) meson is identified by its decays to J/-0(— > e + e )ir + ir ,J/ip(^^i + fi )ir + ir , e + e 
and /i + which have measured branching ratios of 1.95±0.17%, 1.95±0.17%, 0.88±0.13% and 
0.77±0.17% respectively pf]. 



The decay products of the ip{2S) are measured in the CT and are required to have pt > 
120 MeV and 20° < 6 < 165°. Particle identification is provided for the leptons by the LAr 



calorimeter, in which energy deposits typical for muons or electrons are observed, and for muons 
by tracks found in the instrumented iron return yoke. It is demanded that the leptons have 
measured transverse momenta greater than 800 MeV. In the case of the four-particle decays of 
the ip{2S), the lepton pair must have a measured effective mass between 2.8 and 3.4 GeV. No 
particle identification is required for the pions. It is required that no charged particles other 
than the decay products of the ip(2S) are observed in the CT or FT. A measured photon-proton 
centre-of-mass energy, W 7P , in the range 40 — 160 (40 — 120) GeV for the dimuon (dielectron) 
sample is demanded (see section ||) . 

Events in which the scattered positron is observed in the backward electromagnetic calorime- 
ter (BEMC in 1994, SpaCal in 1995) or in the LAr calorimeter, are rejected, which restricts the 
acceptance to the kinematic range Q 2 < 4 (1.2) GeV 2 for the 1994 (1995) data sample. There- 
fore, in section^], positron-proton cross-sections are quoted for Q 2 < 4 GeV 2 . Events with higher 
Q 2 are estimated to contribute less than 1% to the sample. 

The mass distributions obtained are shown in figures ||(a-d), where the data from 1994 and 
1995 and for both di-lepton types are combined. Figure ||(c) shows a clear signal of 27 events with 
negligible background in the range M^s) ~ 60 < AM < Mm^S) ~ MeV where 

AM = M^+^-^+tt-) — Mf£+£-y The number of events in the ip(2S) — > £ + £~ signal is estimated 
by fitting the mass distribution to two gaussian functions, convoluted with an exponential tail 
to account for di-electron decays in which one of the decay electrons has lost energy through 
radiation, above a linear background (figure^(d)). A signal of 50 ± 12 events is obtained. No 
contribution from ^(3770) to these signals is expected. 



5 Monte Carlo Model and Acceptances 

Simulation studies were made to obtain the acceptances and efficiencies of triggering, track 
reconstruction, selection and lepton identification for the tp(2S) production processes. These 



studies used event samples generated with the DIFFVM Monte Carlo generator |14|, and a 



detailed simulation of the HI detector response based on the GEANT program [15]. DIFFVM 
generates events according to the cross-section dependence da/dt oc W^ p e bt for elastic tp(2S) 
production, and d 2 cr/didMy oc W^e^My for production with proton dissociation, where My 
is the effective mass of the hadrons produced in the dissociation of the proton which is related to 
the elasticity z by z ~ (W 2 — My)/(W 2 — m 2 ). The parameters of the DIFFVM generator are 
chosen such that the features of J / ip elastic and proton dissociation photo-production processes 
are reproduced |, g: 4e = 0.9, b = 4.0 GeV" 2 , b' = 1.6 GeV" 2 , = -2.16. Equal cross- 
sections for elastic and proton dissociative photo-production of i/j(2S) are assumed, consistent 
with the results of the J /ip analysis. The decay angular distribution for ip(2S) decaying directly 
into two leptons is simulated according to s-channel helicity conservation. The Monte Carlo 
program is found to model all features of the data well, in particular the momentum and angular 
distributions of the ip(2S) decay products to which the acceptance is sensitive. 

The acceptance is found to be not strongly dependent on W 1P in the range 40 < W 7P < 
120(160) GeV for the electron(muon) decay channel. By demanding that the measured value of 
W 1P lies in this region the sensitivity of the measured cross-section to the uncertainty in its W 1P 
dependence is small. The requirement that there be no charged particle observed in the detector, 
other than the decay products of the ip(2S) meson, restricts the acceptance approximately to 
the region z > 0.95. Events with z < 0.95 are estimated to contribute a background of less than 
1% to the signal. Within the kinematic region z > 0.95 and 40 < W 1P < 120(160) GeV the total 
probability of observing a produced ip{2S) meson is about 6% for the four-particle decays and 
18% for the two-particle decays. 



Uncertainties in the acceptance calculations contribute to systematic uncertainties in the 
cross-section measurements. The uncertainty due to the use of the Monte Carlo model is esti- 
mated to be 5% by observing the changes in the acceptance resulting from varying the parameters 
Ae by ±0.1, b and b' by ±1.0 GeV -2 , (5 between —2.00 and —2.50, and by changing the assumed 
fraction of the quasi-elastic cross-section due to elastically produced ip(2S) by ±20%. Further 
systematic uncertainties arise due to uncertainties in the precision of the simulation of the HI 
detector response. The most significant are 12% due to the track reconstruction efficiency (6% 
for the two-particle decays), 6% due to the lepton identification efficiency and 6% due to the 
trigger efficiency. Additionally for the two-particle decays, there is a 20% systematic uncertainty 
in the cross-section measurement resulting from the choice of the functional form of the fit to 
the di- lepton mass distribution in figure ||(d). The total systematic uncertainties in the cross- 
section measurements are 16% for the four-particle decays and 23% for the two-particle decays. 
Branching ratio (BR) uncertainties of 9% and 18% respectively (l^] are given separately below. 



6 i/j(2S) Cross-sections 

The numbers of events, acceptances, flux factors and cross-sections are given in table |. For 
the four-particle decays of the ip(2S) the results are shown separately for the electron and 
muon decays, and for the 1994 and 1995 data samples. The corresponding four cross-sections 
are consistent with one another. A decomposition of the acceptance is given, which shows 
that, except for the trigger efficiencies, aspects of the acceptance are stable from year to year 
and decay channel to decay channel. The higher trigger efficiencies in 1994 with respect to 
1995 reflect a higher fraction of events in which the trigger fired being recorded for subsequent 
analysis. This fraction is determined by HERA beam background conditions and global data- 
taking requirements. For the two-particle decays all data are combined in table |l] because of the 
lack of statistics of the subsamples in the presence of substantial background. 

The combined cross-section from the four-particle decays using 1994 and 1995 data is 
a [e + p -> e + ip(2S)Y] = [2.02 ± 0.39(stat) ± 0.33(syst) ± 0.18(BR)] nb 
for the kinematic range 40 < W 1P < 160 GeV, Q 2 < 4 GeV 2 , z > 0.95. 

The photo-production (Q 2 = GeV 2 ) cross-section at W 7P = 80 GeV is related to the 
electro-production cross-section by 



(7, 



cp 



a^W B = »G.V)/*/d«J'/ r/ .(,,g»)( 55 ^) l + m 



V>(2S) , 



where W 1P = m 2 — Q 2 +y{s — m 2 — m 2 ), Ae = 0.9, n = 2. The limits of integration are determined 
by the kinematic range of the electro-production cross-section. The spectral flux of transverse 



photons from the electron [16] is 

2m 2 ?/ 2 



/ 7 / e (y,Q 



2\ _. a 



2 vr Q 2 



Q 2 



(3) 



Given an approximate WS p dependence of the photo-production cross-section, the choice of 
W-y P = 80 GeV as the point at which to specify a lp minimises the uncertainty in a lp due to 
possible deviations from this dependence. Varying the assumed W 1P and Q 2 dependences of the 
virtual photon-proton cross-section in equation ^ (4e between 0.8 and 1.0 and n between 1 and 
3) results in a systematic uncertainty in the photo-production cross-section of 3%. 



For z > 0.95 at W 1P = 80 GeV, the photoproduction cross-section 

a [yp ip(2S)Y] = [16.9 ± 3.3(stat) ± 2.7(syst) ± 1.5(BR)] nb 

is obtained. 

The ip(2S) signal of 50 db 12 events from the analysis of the two-particle decay sample yields 
a statistically independent measurement of the cross-sections for ip{2S) production: 

a [e + p etp(2S)Y] = [2.65 db 0.63(stat) db 0.61(syst) ± 0.48(BR)] nb, 

a pyp -» ^(25)F] = [22.1 ± 5.3(stat) ± 5.1(syst) ± 4.0(BR)] nb 

for the kinematic intervals above. Here the estimate of the level of background under the ip(2S) 
signal contributes a 20% systematic uncertainty. 

The results from the samples of two-particle and four-particle decays agree well, so they are 
combined to give an overall photo-production cross-section of 

a [tp ip(2S)Y] = [17.9 ± 2.8(stat) db 2.7(syst) ± 1.4(BR)] nb 

at VF 7P = 80 GeV, for z > 0.95. 



7 The Cross-section Ratio cr [ij>(2S)]/<r [J/i/>] 

The J/ij) quasi-elastic [z > 0.95) photo-production cross-section at W 7P = 80 GeV is (119.5 =b 
11.2 ± 12.0) nb g, yielding a ratio: 

T^^J^lT = 0,150 ± °- 027 ( stat ) ± 0.018(syst) db O.Oll(BR). 

Figure § contains a compilation of measurements from fixed target experiments of the ratio of 
the photon-nucleon cross-sections for the production of ip(2S) and J/ip mesons: three photo- 
production experiments at Q 2 = GeV 2 (deuterium target fl7| SLAC, lithium target (T^] NA14, 
and deuterium target [O] E401) plus muo-production experiments on an iron target pC[|(EMC) 



and more recently on tin and carbon [21|(NMC). Of these the deuterium target data samples 
were restricted to events in which the decay products of the vector mesons were the only tracks 
detected; in the lithium and iron target data some contribution from the kinematic region of 
lower z was included. Despite the different experimental conditions, it is apparent that in 
going from fixed target to HERA energies there is no significant change in the ip(2S) to J/tp 
cross-section ratio. 



8 Colour Dipole Calculation 

By considering the cc quark pair as a colour dipole of transverse size r, cross-sections for the 
photo-production of charmonium can be calculated in the framework of QCD |7], M. The ampli- 
tude of forward photo-production on the proton can be written 

J d 2 r # v (r)V(r)* 7 (r) (4) 

where ^ 7 (r) is the wave- function representing the probability of finding a cc dipole of size r = \ f\ 
in the photon, <j{r) is the colour dipole cross-section, and ^'v'(r) the transverse wave- function of 
the charmonium state. The product a{r)^^{f) is calculated in QCD using the approximation of 



a two-gluon exchange between the colour dipole and the proton. It is found that, approximately, 
x I' 7 (r) oc exp(— m c r) and a(r) ~ r 2 , so that the product is a function peaking at r rs 2/m c . As a 
result the matrix element, and hence the cross-section, is sensitive to the radial wave- function of 
the charmonium around r = 2/m c . This feature influences the ratio of the cross-sections for the 
production of ip(2S) and J/ip(lS) mesons. There is a relative suppression of the ip(2S) because 
its wave- function has a radial node (the value of r for which ^ / y(r) = 0) near to r = 2/m c , in 
contrast to the J ftp wave- function. Therefore the ratio of these cross-sections provides a test of 
the model which is free of many normalisation uncertainties. 

A prediction for the ratio of the forward elastic photo-production cross-sections of 0.17 is 
obtained ||, which agrees well with the ratio of quasi-elastic cross-sections reported here. This 
agreement supports the hypothesis that the charmonium production cross-section is sensitive to 
the radial wave- function of the meson away from the origin, resulting in a suppression of the 
production of the 2S state, ip(2S), with respect to the IS state, J/tp. 

9 Summary 

Measurements of the quasi-elastic photo-production cross-section for ip(2S) at a W 1P value of 80 
GeV at HERA have been reported. The ratio between the quasi-elastic ip(2S) and J/ip photo- 
production cross-sections is compatible with those measured at lower energies. A prediction from 
a QCD based model in which the photon fluctuates into a cc quark pair which then behaves as 
a colour dipole in its interaction with the proton is in good agreement with the measured ratio. 
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Table 1: Summary of signals, acceptances, and the quasi-elastic (z > 0.95) electro- and photo- 
production cross-sections for the five data samples. The first error shown is statistical and 
the second systematic. Uncertainties due to the imprecisely known branching ratios of the 
ip(2S) decays are not included. The ranges over W 1P and Q 2 over which the ep cross section 
measurement is made are shown. Angular and pt track acceptance refer to the fraction of 
events in the kinematic range in which all decay products fall within the angular and pt ranges 
of acceptance, respectively. Trigger efficiency indicates the probability that such an event is 
accepted by the first level trigger, and the analysis efficiency gives the probability of subsequently 
meeting all further selection criteria (see section |J). <I> 7 / e is the photon flux factor used to give 
photo-production cross-sections at a fixed W w of 80 GeV. 
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Figure 3: The effective mass plots for the four track sample: (a) the four track effective mass, (b) 
the di-lepton effective mass, and (c) the AM distribution where AM = Mn+;- 7r + 7r -\ — Mn+i-y 
(d) Effective mass distributions for the two track sample: the di-lepton effective mass with a fit 
to two gaussians, convoluted with an exponential tail to account for di-electron decays in which 
a decay electron has lost energy through radiation before entering the detector, plus a linear 
background. 
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Figure 4: The ratio of the cross-sections for i/j(2S) and J/ip photo-production on nucleons as a 
function of photon-proton centre-of-mass energy, W 7p . The results from the present analysis are 
shown together with those from earlier fixed target experiments [jD], [l^, [l9|, p2l]| . The ratios 
from the fixed target experiments have been corrected using the latest decay branching ratios. 
The error bars shown on the low energy measurements are statistical only. The combined 
statistical and systematic errors are shown for the HI result. Horizontal bars indicate the 
ranges of W~ p over which the measurements were made, where they are specified. The common 
uncertainty due to the imprecisely known branching ratios is not shown. 



